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Abstract 

When charged current weak interations are excluded, the neutral current 
weak interaction is formally similar to ordinary electromagnetism with a 
massive photon. In this spirit, the Maxwell equations for the fields of the 
Z-boson are derived from the standard model. These describe the Z-boson 
scalar and vector potentials, and the Z-boson electric and magnetic fields 
whose sources are electron and neutrino distributions and currents. The Z- 
boson Maxwell equations are solved for point sources representing classical 
point-like electrons and neutrinos. The parity violation of the weak interation 
is manifest in the structure of these solutions. As an application of this 
model, the neutral current contribution to the muonium hyperfine structure 
is computed using nonrelativistic perturbation theory. 
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I. INTRODUCTION 



The electromagnetic force is known to be the result of the exchange of 
photons between charged particles. Outside of the context of quantum field 
theory, this interaction is described in terms of classical scalar and vector 
potentials (A and A) with their electric and magnetic fields (E and B). 
This description of the electromagnetic interaction is a very useful model and 
forms the foundation of a large body of macroscopic, atomic, and condensed 
matter physics. 

The electron also participates in weak interactions, both charged current 
(W-boson exchange) and neutral current (Z-boson exchange). For neutral 
current events, electrons (or neutrinos) remain as electrons (or neutrinos). 
We may imagine that one particle is the source of neutral current electric 
and magnetic fields and that the other particle feels the forces produced by 
these fields, in complete analogy with ordinary electromagnet ism. In the 
charged current case, an initial electron state emerges as a final neutrino 
state or vice-versa. For this reason it is difficult to consider such a picture 
for the charged current interaction. Charged current interactions will not be 
considered in this paper, and hence, we will develop an incomplete picture 
of the weak force. This picture will, however, be useful when an interaction 
cannot proceed by W-boson exchange.. 

This paper begins with an analysis of the standard model Lagrangian for 
the first generation of leptons (the electron and its neutrino). A nonrelativis- 
tic weak-field Hamiltonian for the electron is developed which allows us to 
compute the interaction energy of an electron in the presence of a classical 
Z-boson field. The Maxwell equations for the Z-boson are then developed. 
In the absence of sources, the Maxwell equations are identical to those of or- 
dinary electromagnetism but with a massive photon. The Maxwell equation 
source terms are derived from the interaction energies for both electron and 
neutrino sources. 

The Maxwell equations derived here can be used to describe the (albeit 
small!) Z-boson field generated by macroscopic or atomic-scale distributions 
of electrons. They may also be used to visualize the Z-boson fields surround- 
ing classical point-like electrons and neutrinos. The classical point particle 
solutions provide an interesting visualization of the parity violation in the 
standard model in terms of a vortex-like magnetic field structure oriented 
with the electron's spin. 
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As an application of the formalism developed in this paper, a calculation 
of the neutral current contribution to the hyperfine structure of muonium 
(/i + e~-bound state) is made using nonrelativistic perturbation theory. Other 
applications might be found in nuclear physics. 

II. THE STANDARD MODEL LAGRANGIAN 

We begin by considering the matter terms of the Lagrangian density for 
the standard model (with massless neutrinos) after spontaneous symmetry 
breaking (following the notation of Quigg&): 

^electron = e R (i^d^ + e^A^ - g' sin 9 W ^Z^)e R 

+e L (i7 M <9 M + vfA^ + ^(#cos# w - g' sin 6 w )YZ^)e L 
-me R e L - me L e R , 

(1) 

C-neutrino = V^d^ - -(gCOs9 W + g' SU1 By/^Z^ V . 

Note that the W± boson fields have been set to zero, decoupling electron 
and neutrino in the matter Lagrangian density. We do not wish to con- 
sider charged current interactions. In addition, the neutrino chirality is con- 
strained such that |(1 + 7 5 )i/ = 0. 

The matter bispinors naturally decompose into left and right handed 
chiralities when using the chiral representations of the Dirac matrices: 

'-((-lo- 1 ) (-.0))- « 

In the chiral representation, the matter equations of motion are 



> ( X ( e ^ j = ^electron ( J (3) 

er ■ ir R - $ B -m \ ( e R 

-m -a- ■ 7r L - $ x J \e L 

id V = HneutrinoV = (~tT • 7T„ - $„) V , (4) 
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where 



7Tm 



— a' sin6 w )Z , 



— iV + eA+-((jfcos6 l vK — <7 
-zV + eA—g 1 sin #n/ Z , 
eA° + - (g cos 6>vf - 5-' sin 9 W )Z° , 
eA°-g'sm6 w Z° , 
— iV — -(g cos^V + </ sin#w)Z , 



(5) 



$ n = --(g cos 9 W + g' sin 9 W )Z° , 
(A A) , = ( Z° Z ) 



III. NONRELATIVISTIC-WEAK FIELD APPROX- 
IMATION 

To assess the nonrelativistic limit it will be more convenient to use the 
Dirac representation! for the Dirac Matrices: 



1 
-1 



"" = (In ill _ n M ■ ( 6 ) 
In the chiral representation, the electron Hamiltonian in ([J) becomes 



a 
-er 



/ iff- 



n 



electron 



|<T • (ir R + 7T L ) \ 



\a ■ (ir R + tt l ) 



9^ • (TTr - 7T L J 



(7) 



|(<l> JJ + $ L )-m / 



In the nonrelativistic limit the electron bispinor may be written in the 
form 



-i{m+E NR )t 



(8) 
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where the nonrelativistic energy E^r <C m. We may solve for the "small" 
component spinor x i n terms of the "large" component spinor 0; whence, to 
lowest order in the ratio Emr/ttl: 

Substituting for x i n the equations of motion allows the identification of the 
nonrelativistic Hamiltonian 

n NR = ^Qo-.^ + ttz)-^-^)) (io) 

+ -<7 ■ (tt r - 7T L ) - - + $ L ) . 

Use the identity (<r • A) (<x • B) = A • B+icr- (A x B) and neglect terms 
quadratic in the fields (weak field approximation) to derive 

n NR = --^v 2 + v, (ii) 

2m 

V = --^-V- (eA+ ] - (g cos 9 w -3g' sin e w )Z] 
2m \ 4 / 

(eA + \{gcosO w - 3g'sin6 w )Z ) ■ V 
2m \ 4 / 

+^cr- |vx ^eA + ~(gcos8 w - 3g' siae w )Zjj (12) 

gcos9 w + g'smB w 

8m v ; 

X 1 
— eA° — ~(g cos 6*^ — 3g' sin 9 W )Z° — — (g cos 9 W + g' sin 6^) cr • Z . 

Note that it is understood that the differential operators operate only on the 
object immediately to its right. 



IV. INTERACTION ENERGY 



The interaction energy ( |T2"D is responsible for the Maxwell equation source 
terms. In the nonrelativistic-weak field approximation the interaction energy 
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for an electron source is just the expectation value of V using the nonrela- 
tivistic electron wavefunction 0. After integration by parts we have 



V = (eA+j(gcos6 w -3g'sinO w )Zj ■ (0 f V0 - (V0 f ) 0) d 3 r 

+ i/ ( V>< ( eA + ^(# cos ^- 3 #' sin M z ^ " (<^W) dh 

_. g cos6 w + g >^e w I zQ V0 _ rf3r (13) 

8m J y 

- j (^eA° + ^(g cos 9 W - 3g' sin 9 w )Z°^j f 0d 3 r 

— -((? cos ^vi/ + (?' sin^vi/) j Z- (0^0) d 3 r . 

The nonrelativistic limit makes no sense for a massless neutrino. The 
neutrino source interaction energy may be found directly from TC neut rino given 
by (§: 

V = ~ {gcos6 w +g'sm6 w )J Z°z/W 3 r (14) 
+- (gcosOw + g'sinOyy) J Z- (v^rrv) d 3 r . 

V. MAXWELL EQUATIONS 

The first order Lagrangian for the photon is 

L = J (-E ■ (VA° + A) - B (V x A) + l - (|B| 2 - |E| 2 )^) dh-V . (15) 

From this Lagrangian we derive the relationship between the potentials and 
field strengths and two of the Maxwell equations: 

E = -VA° - A , B = V x A , (16) 

V E = -e0 f , (17) 
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Vx ( B+ ^ - E = - e (2^) - ( v ^) • ( 18 ) 

The first order Lagrangian for the Z-boson field neglecting three and four 
boson couplings (see Appendix A for further discussion) is 

L = J (-E- (vZ° + z) -B-(V x Z) + i (|B| 2 - |E| 2 )^d 3 r 

+ \M 2 j ((Z°) 2 -|Z| 2 )rf 3 r-\/. (19) 

This differs from the photon's Lagrangian by the addition of a mass term 
(M = M z ). Note that the symbols "E" and "B" are now used as the field 
strengths for the Z-boson instead of the photon. The electron source terms 
are modified for the Z-boson: 

V-E+M 2 Z° = -^(gcaa9 w -Zcfwa.0 w )<p<l> (20) 
+^ (g cos 9 W + g' sin 6 W ) (^^j (0 f cr • - ( V0 f ) • acf) , 

Vx ( B+— (gcos9 w -3g' sin W )- <T(f)) ) + M 2 Z — E 
\ 4m 2 ' J 

= - l -( gC os6 w -?>g'sin6 w )(^£^ (0^0-^)0) (21) 

+- (gcos^ + , sin0 W r) - (0 f cr0) . 
Similar Maxwell equations for a massless neutrino are 

V • E+M 2 Z° = - (g cos 9 W + g' sin 9 W ) v^v , (22) 
V x B + M 2 Z-E = -(#cos# w + 

g' sin 9 w ) - (y^crij) . (23) 

The remaining two Maxwell equations are just Bianchi identities and 
remain unchanged from ordinary electromagnetism (Appendix A): 

V • B = , V x E + B = . (24) 
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VI. ELECTRON FIELDS 



It is interesting to consider the fields of a classical point-like electron. 
First consider the photon fields. The wavefunction bilinear combinations 
appearing as sources for the photon electric and magnetic fields have the 
following interpretations: 



— ecjrcj) — > charge density p , 
— e ^2 _ ~) (^^^ ~ i^^) 0) — y current density j , (25) 



cr(f) — ► spin angular momentum density s . 



1 
2 

The term — ^^cf^acf) in (|T8| ) has the physical interpretation of a magnetization 
(or a dipole moment density). From these considerations we find that the 
point-like classical electron has a charge -e and a magnetic moment of — — 
times its spin (a gyromagnetic ratio of 2). 

Now consider the Z-boson fields. There is the addition bilinear form 
(im) {. ( t^ cr ' ^0 ~ i^^) ' <T( / ) ) which may be interpreted as a density for the 
helicity- 1 velocity | product but can be considered vanishing for the electron at 
rest. The Z-boson Maxwell equations for the point-like electron at rest with 
spin polarization (S = ±|) along the z-axis derived from fl2"0|) and (|2T| ) are 
then 

V-E+M 2 Z° = -~(gcos6 w -3g'sm9 w )5(r) , (26) 



Vx (B+—(gcos6 w -3g'smO w )Sz6(r)) + M 2 Z (27) 
\ 4m / 

= - [g cos 6w + 9 singly) Sz5 (r) . 

Unlike the electromagnetic Maxwell equations for the electron at rest, there is 
a nonvanishing magnetic source current density in (p7|) as well as the familiar 
electric charge density in (|26|). 

When the electron is placed in motion along the z-axis in either a state of 
positive or negative helicity, additional source terms are required to maintain 
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the covariance of the Maxwell equations. Specifically, a magnetic source cur- 
rent is generated due to the motion of the electric charge density just as with 
ordinary electromagnetism; however, unlike ordinary electromagnetism, an 
additional electric charge density coupled to the nonvanishing at-rest mag- 
netic source current is required for covariance. Thus, in the limit of a classical 
point-like electron, we must make the identification 

(4^-) (0 f cr • V0 - (V0 f ) • cr0) — ► Sz ■ v5 (r - vt) , (28) 
and the Maxwell equations (|26| ) and ( p7|) become more generally: 



V-E+M 2 Z° = -- (g cos 9 w -3g' sin 6 w )S(r- vt) (29) 
+ - (g cos 6w + g' sin 9 W ) Sz ■ v5 (r — vt) , 



Vx ( B+- — (g cos 9 W - 3g' sin 9 W ) SzS (r - vt) J + M 2 Z 
\ 4m y 

= - (gcos^ + g'sin^) Sz;5 (r — vt) (30) 
— -(gcos#jy — 3g' sin 9 w )v5 (r — vt) . 

The Maxwell equations for the point-like electron at rest (|26| ) and ( |27| ) 
are solved in Appendix B. Physically there are three contributions to the 
Z-boson fields. The first component is in the form of a Yukawa potential due 
to a neutral current electric charge q = — k(gcos9w — 3g' sin 9\y)' 



^Yukawa ^ 



F - q ( M 4- M p~ Mt v 

^Yukawa — — I 1 n 6 V 



(31) 

\ Mr- 

Yukawa . I t „ 

Arc \ r r l j 

The second component is due to a neutral current magnetic dipole mo- 
ment n = — S; in spherical coordinates: 
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W = £(^ + ^)e- Mr sin^, (32) 



H (M 1 



B Dipole = ^^_ + _j e — cos^r (33) 

fi ( M 2 M 1 \ Mr . ~ 2 _ . . 
+f- + — + - e" A/r sin 00 + -/xz5 r • 

Note that (|33| ) implies that the electron has a neutral current gyromagnetic 
ratio identical (when quantum filed theory radiative corrections are ignored) 
to the electromagnetic. 

The third component is a vortex-like magnetic field near the equatorial 
plane due to a magnetic vortex charge kS = \ (g cos 9w + g' sin 9w) S: 



kS 



jpv 



+ ) e- Mr cos#? 



2tt \Mr 2 M 2 r 3 

+ J_ + _JL ) e -Mr gin m ,, , , 

4tt Mr 2 M 2 r 3 ' 

1 kS . 



B PV = ^(M + i_ )e-"'Bin90. (35) 



47T \ r r 



Note that (34) is not finite in the limit M — > 0. This is because there are 



no nonvanishing solutions to B = V x Z and V x B = with the same 
symmetries as (|34]) and fl35l) . The fields fl3"4"D and (|5|) have no direct analogy 
in ordinary electromagnetism, although the axial symmetry of the magnetic 
field ( j35[) is similar to that of a charge in motion. 

The presence of the ^-functions in ( |33"D and (|34|) (discussed further in 
Appendix B) are required topologically to preserve Stokes integral identities. 

The parameter k is a direct measure of the parity violation in the standard 
model since it is the difference between the left-handed and right-handed 
electron coupling constants, while the charge q is their average. Under parity 
inversions — > 7r + 0, 6 ^ n — 6, and r — > :r, 6 — > —0, </>—>• 0. For the 



Yukawa part (31) both Z° and E are even under parity inversion. For the 
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dipole part (|3"2"|) and (|3"3"| ) Z is even and B is odd; however, for fl3~j|) and (P3|) 
Z is odd and B is even, signaling parity violation. 

The sign of the charge q may be freely chosen, using sin 2 9w = 0.23, we 
have |e| = 0.30, \q\ = 0.014, and k = —0.36 sgn(q). 



VII. NEUTRAL CURRENT CONTRIBUTIONS TO 
THE MUONIUM HYPERFINE STRUCTURE 

The /i + e~-bound state (muonium) is an important system for high pre- 
cision studies of quantum electrodynamics (QED).i Its spectrum is qualita- 
tively similar to that of Hydrogen. Using the formalism in this paper, it is a 
simple exercise in nonrelativistic perturbation theory to compute the energy 
level shifts of muonium due to the neutral current interaction. In particular, 
we calculate the neutral current perturbation to the ground state energy of 
muonium, and the hyperfme splitting between the J = (singlet) and J = 1 
(triplet) ground states.^ 

Equation ([13]) is the energy perturbation when (= if) x spinor) is the 
electron's unperturbed nonrelativistic wave function, and Z° and Z are the 
potentials generated by a point-like antimuon. The unperturbed ground state 
electron scalar wavefunction is 

V = ~k=^ , (36) 



V 



na 



where a = 4ir/ (m' e e 2 ) is the muonium Bohr radius with m' e being the reduced 
mass of the electron in muonium. The potentials generated by the antimuon 
are given by equations (|3~TD through (|35D with q replaced by —q, k replaced 
by —k, and the electron's mass replaced by the muon's. Since both particles' 
spins may be freely oriented, it is necessary to replace the fixed choices for the 
polarization axes with Pauli matrices, and in the end, take the expectation 
value using a spin state describing both particles. Following this procedure, 
there are three nonvanishing terms in ([13]): 

V c = q [ Z<Wd 3 r = — ^ ? , (37) 

7 vra 3 (M + 2j of V 1 

jL / ^ R . t . ,/it„ 



V NPV = — — / < B ■ cr > r (38) 

2m ' 
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V PV = --k J < Z ■ <r > (39) 

3vra3 < ^ e '^ > ^( M + 2/a) 2 + M2 

The "cross-terms" in (|i~3D (proportional to Kg) involve the expectation 
value < S e x S M > which can be shown to be vanishing. 

The expectation value of the electron and muon spin operator dot product 
< S e ■ > is evaluated by considering the expectation of the total spin- 
squared S 2 = (S e + S^) 2 , using the decomposition of the symmetric triplet 
spin states, and the antisymmetric singlet spin state into individual particle 
spinors: ||, (|j + IT) jj (triplet), and (|l - ||) (singlet). In the 
triplet state < S e • S M >= |, and in the singlet state < S e • S M >= — |; hence, 
we may compute the neutral current contributions to the hyperfme splitting 
of the ground state of 1.7 Hz from fl3"8|), and 65 Hz from (|39"D, for a total of 
67 Hz. This compares to a 4.5 Ghz electromagnetic contribution which may 
be estimated from ( |38| ) with q replaced by e and M = 0. The contribution 
of ( j37|) shifts both singlet and triplet ground states by -0.10 Hz. 

The neutral current contributions to the spectrum of hydrogen! could 
be analyzed in similar fashion with suitable generalization to include the 
neutral current coupling of the proton. This coupling may be identified by 
summing the contributions from the constituent quarks with inclusion of 
radiative corrections.! Extension to include heavier elements is limited by 
the complexity of atomic wavefunctions. However, for certain cases, such 
as Cesium, accurate calculations can be performed.^ Also of interest are 
muomc aton J because the muon's Bohr radius is much smaller than the 
electron's. In all these systems, the neutral current interaction introduces 
parity violation into the atomic structure due to the mixing of parity-even 
and parity-odd pure electronic states. 

The analysis of positronium (e + e~-bound state) is , however, more prob- 
lematic because of the annihilation orocess e + e~ — > virtual photon — > e + e~ . 
A precise treatment of positroniuma requires using the Bethe-Salpeter equa- 
tion for relativistic bound states. 
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VIII. NEUTRINO FIELDS 



The Maxwell equations for a point-like massless neutrino travelling through 
the origin at t = towards the positive z-axis, after making the identifications 



— > 5(r — zt) and ~ (z/^cri/) — > Sz5(r — zt), are 



wis 



V ■ E+M 2 Z° = n5(r - zt) , (40) 



V x B + M 2 Z - E = -2/«5z5(r - zt) . (41) 

Appendix C considers the solution of the Maxwell equations (fiOD and 
PH). A physical (5 = —5) point-like massless neutrino travelling along the z- 
axis at the speed of light generates scalar and vector potentials nonvanishing 
only in the plane perpendicular to z at the neutrino's instantaneous position. 
An observer could then only detect the neutrino's presence at its moment 
of closest approach. This is not a property unique to the neutral current 
force; rather, it is a property related to Lorentz transformations properties 
of potentials and field strengths. The same phenomenon may be found in the 
ultra-relativistic limit of the Lienard-Wiechart potential for the electron.!^ 
The neutral current fields are given in cylindrical coordinates below: 

Z° = £-K (Mp)6(z-t) , Z = ^K (Mp)5(z-t)z , 

Z7T Z7T 

(42) 

E = ^L Kl (Mp)6(z-t)p, B = ^-K 1 (Mp)S(z-t)4>. 

Z7T Z7T 

where k is the neutrino's neutral current charge, and Kq and K\ are modified 
Bessel functions. 

Under parity inversion the fields ( pffi) represent a solution propagating 
backwards in time (antiparticle). Under both parity and time-reversal, E 
and B are even, while Z° and Z are odd. 

If we consider solutions to the equations (|4"0"D and (ffl]) for neutrinos of 
the wrong helicity (S = +|), the formal solution to ( fiOf) and ( ffl[) would then 
contain delta function derivatives (Appendix C). 
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IX. CONCLUSIONS 



A nonrelativistic weak-field Hamiltonian suitable for describing neutral 
current interactions has been derived from the standard model. Maxwell 
equations for the Z-boson fields have also been derived. They are suitable 
for a description of the small Z-boson fields at an atomic or macroscopic 
scale. 

The fields surrounding classical point-like leptons due to the neutral cur- 
rent interaction have been derived from the Z-boson Maxwell equations. The 
electron possesses Z-boson electric and magnetic fields analogous to its elec- 
tromagnetic Coulomb and magnetic dipole fields. It also exhibits a Z-boson 
magnetic field forming a vortex-like structure oriented with the spin. This 
field is an explicit manifestation of the parity violation of the standard model. 
The fields of a point-like neutrino have also been shown to be nonvanishing 
only in the plane normal to its spin and containing the particle. 

Using the nonrelativistic Hamiltonian and the solution of the Maxwell 
equations, we have calculated the 67 Hz contribution to the hyperfine split- 
ting of the muonium ground state using simple nonrelativistic perturbation 
theory. 
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XL APPENDIX A: FIELD EQUATIONS FOR THE 
STANDARD MODEL 

The gauge terms in the standard model Lagrangian are 

C - — I f fv v 

^ gauge — ^ a^v 1 a ^J/J-vJ i 

Fa^v = d v b a ^ — d^bau + g£ a b c bb^b cv , (Al) 
f = 8 A - 8 A 

where A^ and b afl are respectively the U(l) and SU(2) gauge fields, linear 
combinations of which form the physical boson fields: 
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Z = -Asm 9 w + b 3 cos O w , (A2) 

A — A cos 9 W + 6 3 sin 9 W . 

Using the definitions (A2), field strengths for the physical fields may be 
defined as 

F% = d v W± - d,W± ± ig cos 9 W (W±Z U - ZJV±) ± xe (W±A U - A^W±) , 

F% = d„Z, - d,Z u - ig cos 9 W (W+WJ - W~W+) , (A3) 
= d u A, - d,A u - te (W+W- - W~W+) . 

where e = g sin 9w = 9' cos Ow- The gauge Lagrangian density now including 
the boson mass terms, becomes 



C„ = -\f%F a ^ - -F^F 2 ^ - ±F+F-r + M W W+W~» + \m\Z^ . 

(A4) 

The equations of motion are 



cfF% ± ig cos 9 W F%Z V ± ieF^ v A v T ieF*W ±v =F ig cos 6 W F^W ±V + M 2 W W+ 
= weak charged current source (e.g., e — > z/ and 1/ — > e) , 

0"F£ + i<? cos W F" - ig cos W F+ W"" + M 2 Z Z, (A5) 
= weak neutral current source (e.g., e — > e and 1/ — > 1/) , 

+ ieF" - ieJ£ 
= electromagnetic source (e.g., e — > e) . 

where the matter terms in the full Lagrangian produce the sources for the 
boson fields in equations (A5). 

To characterize the classical fields surrounding a lepton we consider only 
configurations for which the lepton remains invariant in its initial state. 
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There are no weak charged currents and hence W ± = is a consistent solu- 
tion to the fields equations and the remaining fields equations for A and Z 
decouple. 

There are also Bianchi identities for the fields strengths which are satisfied 
identically: 

cT *F± ±ig cos0 w * F%Z» ±ie* F±A v ^ie *F*W ± ^ig cos 9 W *F Z W ± » = , 



d u *F^ + ig cos 6 W *F~ V W +V - ig cos 6 W *F+W~ U = , ( A6) 

gr *f£ + ie *f- u w + " - ie *f; v w~» = , 

where *F IXV = ^e^ upa F pa . When W ± = the Bianchi identities decouple, 
and since there are no boson mass terms in the Bianchi identities, they reduce 
to the familiar form of classical electrodynamics for both the photon and the 
Z-boson. 

The structure of the classical field theories for A and Z are now identical 
except that the Z-boson has a mass M = Mz which appears in its equation 
of motion but not in its Bianchi identity. 



XII. APPENDIX B: SOLUTION FOR THE ELEC- 
TRON 

Without sources, B and Z are solutions to Vx (V x U) + M 2 U = 0. 
When U is a solution then so is V x U. Two solution pairs in spherical 
coordinates are 

U =*"M J , \ sinfl</>, (Bl) 

3 r is (z) i ^ 3 r zh ( z ) - is ( z ) ) 

\Ks_{z)\ \-zKi{z)-Kz(z)\ 

where I n (z) and K n {z) are modified Bessel functions with z = Mr; 
explicitly 
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Il(z) = \ — ( coshz - | , Ii(z) = \ — sinhz , 

2 W ]/ nz \ z J 2 V nz 

K i^ = \fl{ 1 + i) e "- K ^ = irr- (B2 » 

Consider a spherical volume V containing the uniform magnetic source 
density ySz — yS ^cos# r — sin# O^j. Consider homogeneous solutions for 
Z and B inside (upper) and outside (lower) the sphere: 

3 ( als (z) ~] _ 3 f azli (z) — alz (z) 1 ~ 



bKz (z) J I -bzKi (z) - bK % (z 

i f ali (z) . 

B = -Mz-2\ *)'}smO<i>. (B3) 
1 bKz (z. 



2 



A particular solution to V x B + M Z =ySz must be added inside the 
sphere, such as 

Z = M~ 2 ^Sz , B = . (B4) 

The constants a and b are determined from the boundary conditions that 
the components of Z and B parallel to the interface at r — R are continuous 
across that interface: 

a = -£-SM (MR)~^ Kz (MR) -> R~ 3 as R -> , (B5) 
6 = -— (Mi?)" 3 / 3 (Mi?) -> -—SM\ - as i? -> . 

47T 2 47T V 7T 

The point-like electron solution corresponds to the R — > limit; however, 
note that the constant a diverges like R~ 3 implying that Z may possess a 
5-function at the origin. To determine its weight, integrate V x B + M 2 Z 
over a spherical volume and apply Stokes' theorem for volume integration: 



kSz = / (V x B)d 3 r+M 2 / Zdh = / (d 2 S x B)+M 2 / Zdh . (B6) 
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This holds for any radius but in the limit as the radius tends to oo the 
magnetic field vanishes. Z must have a term (r)z for equality to hold. 

Next consider a spherical volume V which contains the uniform magne- 
tization M =-(jZ = ^ I cos 6 r — sin 6 6). Solutions for Z and B are 



bKs(z) 



sin0<£, (B6) 



; I ali (z) 1 _ 3 f azh (z) — a/3 (2) 1 

B = -2Mz~^i * y 1 \cosdv + Mz'i I * v ' * K ' \sixi60 
Mi{z)\ \-bzKi_(z)-bKz{z). 



2 2 



The vector potential Z must be continuous across the interface at r = R. The 
other boundary condition involves H = B — M. It is easy to demonstrate 
using Stokes' theorem that, without surface current sources, the tangential 
component of H must be continuous across the interface, The solutions for 
the constants a and b are 

a = ~^M 2 (MR)^ Ks (MR) -> iT 3 as R -> (B7) 



6 = (M#)~5 /a (MR) -> y - as i? ^ 

The divergence of the constant a signifies the presence of a 5-function in 
the magnetic field. To determine its weight, compute the volume integral 
of V x Z = B and use the Stokes theorem for volume integration as be- 
fore. There is a similar singularity in the point-like electron's electromagnetic 
dipole field.0 

The solution for the electric field is a straightforward application of Gauss' 
law. 



XIII. APPENDIX C - SOLUTION FOR THE MASS- 
LESS NEUTRINO 

Consider a massless neutrino whose trajectory takes it through the origin 
and then along the positive z-axis. Its potentials and fields will be functions 
of u = z— t and the cylindrical coordinate p. The Maxwell equations fl40|) and 



(|4T|) may be written as second order differential equations for the potentials 



using cylindrical coordinates: 
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B2Z+ ldZ \ M ^-2^ + im + l -Z>) = £(l-2S)«(r-*) 



<9p 2 p dp du 2 du \ dp p J 2 

d 2 Z- ldZ~ .,„„ re 



(CI) 



+ M 2 Z~ = - (1 + 25) 5(r - zt) (C2) 



<9p 2 p <9p 2 

S-5 Z ' = ° (C3) 



dp 2 p dp \p 2 

where Z + = \ (Z° + z-Z), Z~ = § (Z° - z-Z), Z' 5 = p • Z, and Z* = </> • Z. 

Equation (C4) is the modified Bessel equation. With no singularity at 
the origin, equation (C4) has no bounded solution at infinity other than 
Z^ = 0. Similarly, for S = -\ we must choose Z~ = Z p = 0. Equation (CI) 
is then satisfied by the modified Bessel function Z + (p,u) = aK (Mp)S(u). 
To determine the weight a, Gauss's law is applied to a tube of radius R 
concentric with the z-axis: 



k = J E ■ d 2 S+M 2 J Z°d 3 r = 27raMRK 1 {MR) + 2iraM 2 J K {Mp)pdp 
= 2na]imMpK 1 (Mp) = 2rra (C5) 

When S = +|, equation (C2) has a solution proportional to Ko(Mp)S(u). 
The 5-function is then differentiated in (CI) and (C3). 
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